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Abstract. The imager on board INTEGRAL (IBIS) presently provides the most detailed sky images ever obtained 
at energies above 30 keV. The telescope is based on a coded aperture imaging system which allows to obtain sky 
images in a large field of view (29° x 29°) with an angular resolution of 12'. The System Point Spread Function 
of the telescope and its detailed characteristics are here described along with the specific analysis algorithms 
used to derive the accurate point-like source locations. The derived location accuracy is studied using the first 
in-flight calibration data on strong sources for the IBIS/ISGRI system. The dependence of the calibrated location 
accuracy with the signal to noise ratio of the sources is presented. These preliminary studies demonstrate that 
the IBIS/ISGRI telescope and the standard scientific analysis software allow source localizations with accuracy 
at 90% confidence level better than 1' for sources with signal to noise ratios > 30 over the whole field of view, in 
agreement with the expected performances of the instrument. 
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1. Introduction 

The IBIS telescope (Imager on Board of the INTEGRAL 
Satellite) (Ubertini et al. 2003), launched onboard the ma- 
jor ESA gamma-ray space mission INTEGRAL (Winkler 
et al. 2003) on October 2002, is a hard X-ray / soft 7- 
ray telescope based on a coded aperture imaging system 
(Goldwurm et al. 2001). The IBIS imaging system and the 
IBIS scientific data analysis are described in Goldwurm et 
al. (2003). Here we discuss more in detail the specific im- 
age data analysis procedures used to evaluate the point 
source locations and we present the characteristics of the 
System Point Spread Function obtained with the imple- 
mented analysis. We also provide preliminary results on 
the source location accuracy obtained from the first in- 
flight calibration data of ISGRI, the low energy (15-1000 
keV) detector of IBIS IjLebrun et al. 2003|l . 
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2. The IBIS image decoding and the associated 
System Point Spread Function 

In coded aperture telescopes (Dicke 1968, Fenimore & 
Cannon 1978, Caroh et al. 1987, Goldwurm 1995), the 
source radiation is spatially modulated by a mask of 
opaque and transparent elements before being recorded 
on a position sensitive detector. Reconstruction of the sky 
image is generally based on a correlation procedure be- 
tween the recorded image and a decoding array derived 
from the mask pattern. For the IBIS system, where the 
mask is built by the replication of a Modified Uniformly 
Redundant Array (MURA) basic pattern (Gottesman & 
Fenimore 1989) of the size of the detection plane, the raw 
image projected on the detector by a source in the Fully 
Coded Field of View (FCFOV) (sky region from where 
the recorded source radiation is fully modulated by the 
mask) will be a shifted version of the mask basic pattern. 
A source in the Partially Coded FOV (PCFOV) wiU in- 
stead project only a part of the mask pattern. The total 
recorded image (the shadowgram) is therefore the sum 
of the shadows projected by all sources in the FOV plus 
a background term. For those masks patterns for which 
a correlation inverse matrix exists, the case of the IBIS 
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Fig. 1. The ISGRI SPSF compared to the analytical approximation used to fit source positions in the FCFOV (left) 

and in the PCFOV (right). The reconstructed source profiles (solid line) before normalization to the FCFOV for the 
Crab and the fitted bidimensional Gaussians (dotted line) are shown (fitted FWHM are 2.69 pixels for the FCFOV 
and 3.41 pix for the PCFOV). Also shown is the source model profile (broken line) computed at the position derived 
from the Gaussian fit. 



system (Goldwurm et al. 2003), the sky image is recon- 
structed by correlation of the shadowgram with a de- 
coding array obtained from its correlation inverse matrix 
(Fenimore & Cannon 1981). For those systems (and with a 
perfect detector) the resulting sky image of a single point- 
like source in the FCFOV will have a main peak at the 
source position, flat side lobes in the FCFOV and cod- 
ing noise with 8 main source ghosts in the PCFOV. For 
a PCFOV source, a main peak will still be present at the 
source position but coding noise can extend all over the 
image starting from the close proximity of the main peak 
and ghosts will appear in the FOV. The distribution of 
the coding noise does depend on the mask pattern used. 
For example URA masks built from quadratic residues 
(Fenimore & Cannon 1978) have a high degree of sym- 
metry along the axis of the mask, unlike those built from 
Hadamard arrays (Proctor et al. 1979). This leads the 
coding noise to concentrate along the image axis pass- 
ing through the source main peak, producing positive and 
negative sidelobes (see Fig. 3 in Goldwurm et al. 2003). 

We refer to the System Point Spread Function (SPSF) 
as the spatial response of the system to a point-like source 
after image deconvolution, i.e. considering the decoding 
process (Fenimore & Cannon (1978). We will discuss here 
the SPSF relative to the type of deconvolution described 
in Goldwurm et al. (2003), used in the IBIS standard 
data analysis and optimized to obtain the best signal to 
noise (S /N) for point-like sources. This is basically a kind 
of balanced finely sampled cross correlation (Fenimore & 
Cannon 1981) extended to the PCFOV (Goldwurm 2001, 
Goldwurm et al. 2003) and in which the sampling of the 
decoding array is performed in such a way to weigh the 
detector pixels with the fraction of transparent or opaque 
area projected by the mask elements (sect. 3). 

The algorithm implemented is fully described in sect. 2 
of Goldwurm et al. (2003). From the mask array M (= 
for opaque elements and = 1 for transparent elements) 



two decoding array G+ and G~ are obtained projecting 
the arrays M and 1 — M over a detector pixel grid and 
padding the zones outside the mask with Os. Then the sky 
image S is derived from the detector image D by appling 
the following operation for each i,j sky pixel and where 
sums run over all k,l detector pixels 

kl kl 

where Ba = ^fc' ^i+k.^+i^''' -g ^j^g balance array and the 

array W is set to for dead or noisy pixels and to 1 for 
active good pixels. The correlation performed is balanced 
in the sense that in absence of sources and with a con- 
stant background term, the reconstructed images are flat. 
The decoded images are then properly renormalized to the 
reconstructed counts in the FCFOV. 

The above deconvolution for a perfect URA coded 
mask system would provide for a source in the FCFOV a 
SPSF close to a pyramidal function with totally flat side- 
lobes in the FCFOV (Fenimore & Cannon 1978). However 
since the coding is not perfect, due to several instrumental 
effects like those produced by dead zones between pixels 
and detector modules and by the supporting structures, 
the reconstructed image will contain source sidelobes. In 
the PCFOV sidelobes are inherent to the decoding pro- 
cess even for a perfect detector system and the recon- 
structed source peak for this decoding procedure can be 
larger or distorted because of the imperfect coding (sect. 
4). Because of the presence of the ghosts and of the coding 
noise, the general sky image reconstruction process is an 
iterative procedure that combine the initial decoding pro- 
cess to the analysis, modelling and cleaning of sidelobes for 
each source detected in the field. In crowded fields the ini- 
tial source localization may be disturbed by the presence 
of other source sidelobes, but as the modelling of the sky 
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is improved during the iterations the locahzation can im- 
prove. However the performances of the locahzation in this 
case depends on the number and distribution of sources in 
the FOV and the performances of the iterative procedure 
including the convergence criteria used. We present here 
the characteristics of the SPSF at the first step of the it- 
erative procedure after the first decoding and before the 
sidelobe cleaning, and we focus on the characteristics of 
the main peak of the SPSF around the source position. 

The main peak angular width will be approximately 
of the projected angular size of the mask element, i.e. for 
the IBIS telescope w 12' (FWHM). However the exact 
shape of the SPSF depends also on the detector spatial 
resolution, the real features of the telescope (e.g. presence 
of dead and off-pixel zones, mask thickness, etc.) and on 
the decoding algorithm. 

3. Analytical approximation and fitting procedure 
of the SPSF 

The pixels of the IBIS detectors are smaller than the mask 
elements. The ratio R of the linear mask element size to 
the linear pixel size is 2.43 for ISGRI. In order to bet- 
ter sample the shadowgram, the corresponding decoding 
array is also sampled at the pixel scale by redistributing 
its values according to the fraction of pixel area covered 
by a given projected mask element. The resulting decoded 
image has pixels of the angular size of the detector pixels 
(« 4.94'). This way of sampling the decoding array before 
deconvolution optimizes the signal to noise ratio of point- 
like sources in the reconstructed image, as it takes into 
account the blurring induced by the finite spatial resolu- 
tion of the detector (discrete pixels) (Cook et al. 1984). 

For this reconstruction the theoretical SPSF peak in 
the FCFOV is space-invariant and given by the convo- 
lution of two square-pyramidal functions with FWHM's 
equal to Wm, mask element size, and Wp, linear pixel 
size, respectively. Using the central limit theorem we see 
that the convolution of these functions can be approxi- 
mated by a bi-dimensional Gaussian with a width of Wgpsf 

« ^ w'^ + Wp . For IBIS/ISGRI, using pixel units, Wp = l, 

Wjn=2A3, Wgpsf — 2.62. This is slightly worse than the 
theoretical angular resolution of one projected mask el- 
ement Win- The standard imaging analysis procedure of 
the IBIS data performs, for each detected significant ex- 
cess in the deconvolved image, a fit between an im- 
age sector around the source peak and a bi-dimensional 
Gaussian. The free parameters are the centroid position of 
the Gaussian (2 parameters), the 2 (variable in PCFOV) 
widths along the 2 axis, the amplitude of the Gaussian 
and a constant level (background). 

In Fig. n we show the reconstructed peak of a strong 
point-like source for the IBIS/ISGRI system, versus the 
best fitted bi-dimensional Gaussian, in both the FCFOV 
(left) and the PCFOV (right). The source is the well 
known Crab (nebula and pulsar) which was observed 
with INTEGRAL in February 2003. The best fit Gaussian 



width is about 2.65 pixels which is compatible with Wspsf ■ 
A bi-dimensional Gaussian function is a reasonable ap- 
proximation of the SPSF peak except for the negative 
wings around the source, due to the non perfect coding 
of the detector plane (in the FCFOV mainly due to the 
dead zones). These side lobes will be corrected using the 
source model once its fine position is determined. Fig. ^ 
shows also the shape of the peak obtained from the de- 
convolution of the source model computed a posteriori for 
the position obtained with the Gaussian fit. One can see 
that the SPSF shape (including the wings) is well repro- 
duced by the model. However, the computation of such a 
model and its deconvolution is time consuming and thus 
difficult to use for fine determination of the source po- 
sition. The source location is therefore determined by a 
fitting procedure with a Gaussian function and the errors 
by the standard computation involving the curvature ma- 
trix (Press et al. 1996). The error computation assumes 
gaussian distribution of the counts and also independence 
between pixels. In such decoded image the sky pixels are 
instead highly correlated in particular on length scales 
of the mask elements. Moreover the residual background 
structures may make the distribution highly non-gaussian. 
While the optimization procedure is still valid, the derived 
goodness of the fit and error determination may suffer 
of these conditions. The computed formal statistical er- 
ror can be underestimated and systematic errors may be 
dominant, therefore it is important to evaluate the uncer- 
tainties using in-flight data and studying the systematic 
effects in the procedure. 

4. The shape of the SPSF in the PCFOV 

Outside the FCFOV the SPSF is not space invariant even 
if the detector plane were perfect, because the optimum 
properties of the MURA mask are not respected. The 
shape of the SPSF we obtain for the decoding process 
applied can be distorted or enlarged (Fig. ^ right). We 
have studied the properties of the SPSF in the total FOV 
of the telescope and the performances and limitations of 
the Gaussian fitting procedure using the source model. 
Fig.Elshows the fitted width of the SPSF on both axis for 
source positions all along the FOV. The width is approxi- 
mately constant in the FCFOV and ranges around 2.6. In 
the PCFOV however it increases up to values of about 4 
pixels (« 20'). The behavior is symmetric with respect to 
the image center and similar along the two axis. As a con- 
sequence off-axis sources can appear elongated or spread 
in the reconstructed images. 

By studying the offset between the fitted and input 
source model position, we have derived a measure of the 
systematic bias introduced by the approximation made for 
the form of the SPSF. The effect is not large (< 0.5') and 
can be estimated. An automatic correction of the bias has 
been included in the localization procedure. From the first 
analysis of in flight IBIS/ISGRI data, a systematic offset 
between the IBIS telescope axis and the axis of the satellite 
star sensors used to reconstruct the absolute spacecraft at- 
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Fig. 2. The variation of the width of the SPSF (for the 
adopted decoding and before sidelobe cleaning) along the 
FOV of the IBIS/ISGRI telescope, for the two axis (solid 
line for the Y axis and broken line for the Z axis). The 
two horizontal lines indicate the widths Wm and Wspsf and 
the two dotted vertical lines the limits of the FCFOV. 

titude was found. This predominant (« 8'-10') effect, was 
measured using the localization procedure described here 
and the available INTEGRAL data. A matrix for mis- 
alignement correction was computed (Walter et al. 2003) 
and included in the localization algorithms which oper- 
ate the conversion between sky pixels and celestial coordi- 
nates. The imaging software has already provided good re- 
sults for the source absolute positioning with IBIS/ISGRI 
data (see this volume for a number of results). 

In the following we characterize the performance of 
the system (the telescope and the software) in term of the 
point source location accuracy obtained with the most re- 
cent version of the IBIS /ISGRI specific analysis software 
that will be implemented in the October 2003 release of 
the ISDC Off-line Scientific Analysis (OSA 3.0). The loca- 
tion algorithms have been improved and provide now lo- 
cation accuracies for off-axis sources comparable to those 
of on-axis sources, therefore with better results than pre- 
viously reported (Walter et al. 2003). 

5. The IBIS point source location error 

From the formal errors computed through the curvature 
matrix, it can be shown that the average Point Source 
Location (statistical) Error (PSLE) for an optimum coded 
aperture system with a defined SPSF depends on the 
source signal to noise ratio (S/N) as PSLE ~ x \s/n) ■ 
The S/N is the ratio between the reconstructed source 
peak and the computed statistical standard deviation as- 
sociated to the peak pixel. Using computation of such er- 
ror and simulations performed on a perfect system with the 
geometrical characteristics of the IBIS/ISGRI telescope, 
we derived the expected PSLE dependence on the S/N 
(Goldwurm et al. 2001). This curve represents the best 
theoretical accuracy of point-like source location which 
can be attained by such a system. 



With the described analysis software we have derived 
the locations of several strong known sources, for differ- 
ent pointings and energy bands. This has provided us 
with a number of measured source locations in a wide 
range of S/N ratios and angles from the telescope axis. 
We have computed the offsets between the derived and 
catalogue source positions and compared their 90% con- 
fidence level dispersion to the theoretical PSLE values. 
Fig-EIleft) shows the results obtained combining all mea- 
sures performed for source S/N > 6. The S/N bin widths 
were defined to include a constant number of measured 
offset values, in order to have similar precision in each 
bin. Sources at any distance from the telescope axis were 
measured and no dependence with axis angle or energy is 
observed. It can be seen that measured offsets are typi- 
cally comprised between 3' and 20" and are better than 1' 
for S/N > 30. Although the points do not exactly follow 
the theoretical PSLE curve, there is a clear trend compat- 
ible with the expected dependence on S/N (see derived 
parameters in Fig. . At high S /N the dispersion reaches 
a constant level of about 20" which shows the maximal 
accuracy obtained. Residual systematic effects (for exam- 
ple due to background structures) may still influence the 
dispersion, and work is in progress to fully evaluate their 
impact on the source location determination. 

In Fig. O (right) we report the error radius versus a 
measured S/N rather than the pure statistical S/N. The 
measured S /N is obtained using a standard deviation com- 
puted directly from the reconstructed sky image. In our 
analysis we measured the standard deviation in a region 
of 2.5° X 2.5° around the source after cleaning the side 
lobes of all sources detected in the FOV and neglecting the 
source peak regions. This way to evaluate the S /N includes 
an estimate of the residual systematic errors. The derived 
offsets are now much closer to the theoretical PSLE curve 
and indicate that an even better evaluation of location er- 
ror may be obtained when systematic effects will be fully 
accounted for. 

We stress that these results involve fields with few de- 
tected sources. For crowded fields the location estimate 
may depend on the efficiency of the iterative cleaning al- 
gorithm and a full detailed evaluation is not yet available. 
On the other hand, we have verified that the location ac- 
curacy obtained by applying the same algorithm and ap- 
proximation of the SPSF to mosaicked images (Goldwurm 
et al. 2003) is further improved. The results we have ob- 
tained for selected sources in combined images show that 
their localization indeed improves for a given statistical 
S/N. This may be due to the fact that the systematic ef- 
fects present in single pointings are smoothed when images 
are summed. Using Cyg X-1 data, offsets of the order of 
5"-10" have been obtained for S/N values > 300. 

These results show that the IBIS/ISGRI telescope cou- 
pled to the analysis procedures we have developed pro- 
vide point-like source locations with accuracies which fully 
comply with the expected performance of the instrument. 



Gros et al.: IBIS SPSF 5 




50 100 150 200 250 300 20 40 60 80 100 1 20 

SNR meosured SNR 

Fig. 3. The 90 % confidence level (c.l.) error radius on the point-like source location of the IBIS/ISGRI telescope vs. 
the source signal to noise ratio (S/N), as derived from in-flight calibrations. S/N bins include 30 measures of the offset. 
The values are compared to the theoretical 90% c.l. PSLE (solid line). The plotted data were fitted with a function of 
the type y — ax'^ -\- b (broken line). Left: values derived from ^ 2000 offset measures of reconstructed Crab, Cyg X-1 
and Cyg X-3 locations in energy bands between 20 and 300 keV, with source positions ranging between 0° and 14° 
from the telescope axis. The best fit parameters are a = 22.1, b = 0.16 and c — —0.95. Right: Same values plotted 
versus a measured S/N (see text). 
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